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Abstract 

We propose a renormalizable dark matter model in which a fermionic dark mat¬ 
ter (DM) candidate communicates with the standard model particles through two 
distinct portals: Higgs and vector portals. The dark sector is charged under a U{iy 
gauge symmetry while the standard model has a leptophobic interaction with the dark 
vector boson. The leading contribution of DM-nucleon elastic scattering cross sec¬ 
tion begins at one-loop level. The model meets all the constraints imposed by direct 
detection experiments provided by LUX and XENONIOO, observed relic abundance 
according to WMAP and Planck, and the invisible Higgs decay width measured at 
the LHC. It turns out that the dark matter mass in the viable parameter space can 
take values from a few GeV up to 1 TeV. This is a new feature which is absent in the 
models with only one portal. In addition, we can find in the constrained regions of 
the parameter space a DM mass of ~ 34 GeV annihilating into h quark pair, which 
explains the Fermi-LAT gamma-ray excess. 
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1 Introduction 


Cosmological observations indicate unequivocally that the Standard Model (SM) 
particles constitute only 5% of the mass content of our Universe, the rest 26% dark 
matter (DM) and 69% dark energy are yet unknown [HIS]- If the weakly interacting 
massive particle (WIMP) is a correct scenario for DM, a big question however, is 
what the dark matter is made of and what would be the fundamental interaction of 
its constituents with the ordinary matter. Direct detection experiments are designed 
to probe dark matter (DM) elastic scattering off nuclei. In this regards, underground 
LUX [3] and XENONIOO [4] dark matter experiments so far have found no signal on 
these type of interactions, even though they provide us with an upper limit on the 
elastic scattering cross section. 

Recent observation of the galactic center gamma-ray excess (GCE), given its 
intensity and spatial morphology, can be explained by the dark matter annihilation 
in the galactic center (GC) [Hbllj. There is a large number of models in the literature 
suggested in order to explain the gamma-ray excess but among them are models with 
dark matter, annihilating predominantly into b quark pair and on top of that can 
evade stringent bounds from direct detection experiments UMi. 

Moreover, there are models as extensions to the SM with an additional U{1)' gauge 
symmetry and its associated Z' boson which contain a DM candidate annihilating 
via an intermediate neutral gauge boson Z' [T9II29] . 

Models with an extra broken U{1)' gauge symmetry are motivated by new physics 
beyond the SM, as examples for various models we recall, those focusing on the grand 
unified theories like 50(10) and Eq, (see e.g. [30l|3T] and for a review consult [32]), 
dynamical symmetry breaking models like topcolor (see |33j for a review), decoupled 
models like leptophobic Z' |3l|, Little Higgs theories [35H38] . Twin Higgs model [39|, 
family non-universal scenario @01. 

There are many models where the SM and DM sectors interact through only one 
portal, in the sense that there is only one type of mediator (e.g. scalar or vector) to 
connect the two sectors. There are also models with two similar mediators @I1|42|. 

In this work we propose for the first time a minimal dark matter model with two 
distinct portals in tandem: Higgs and vector ones. In this article we construct a 
two-portal DM model based on an extra gauge symmetry 1/(1)' which not only can 
explain the observed relic density and the galactic gamma-ray excess but also can 
evade direct detection as well as the constraints from invisible Higgs decay. 

In our two-portal model the dark sector consists of a Dirac fermion as a WIMP 
dark matter candidate, a complex scalar field as the first mediator, both of them 
charged under a new 11(1)' gauge symmetry. Obviously the fermion dark matter is 
coupled to the dark gauge boson Z', covariantly. In addition, we assume that only 
the SM quarks are also charged under the 17(1)'. Thus, the Z' field interacts with the 
SM particles and the DM, hence the second mediator of the model. In the present 
model we will deal with a non-universal Z' gauge boson, such that the new gauge 
boson has negligible coupling to the first and second generation of the SM quarks. 

Models with a new Z' gauge boson which prefer interaction with only the third 
family of the SM fermions are widely investigated within different scenarios beyond 
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the SM. Among these scenarios, we recall superstring inspired models [33], topcolor 
assisted technicolor model |43j . phenomenology of a Z' boson coupled only to third- 
family fermions [44j . electroweak constraints on models with non-universal Z' bosons 
|45j . and warped models in which the extra Z' boson typically couples only to the 
third generation |23) . 

The present article has the following structure. In the next section we introduce 
our model in detail. In Sec. [3] we compute the SM Higgs invisible decay width within 
the model and address the constraints on the invisible decay width from the LHC 
measurements. We derive a formula for the DM-nucleon cross section in Sec. m In 
Sec. [5] we discuss on the DM relic density within the thermal freeze-out mechanism 
and our numerical computations for the relic abundance and DM-nucleon interaction 
are discussed in Sec. m In Sec. El we will find regions in the viable parameter space 
of the two-portal model that can explain the Galactic gamma-ray excess given the 
recent Fermi-LAT data analysis. We finally finish with a conclusion. 


2 The Model 


We introduce a dark matter model which has the property of having DM-SM inter¬ 
action through two different portals, i.e. the vector portal and the Higgs portal. We 
will show in Sec. Hjthat for the model to be elaborated below there is no tree level DM 
scattering off nuclei and the Feynman diagrams begin with a one-loop contribution 
which turns out to be suppressive. A two-portal dark matter model can therefore 
be designed in order to predict the DM elastic scattering to be consistent with the 
direct detection experiments. This is reminiscent to the velocity suppressed models 
for elastic scattering processes but in a different way. 

The details of the model comes in the following: beside having a scalar field that 
mixes with the SM Higgs via the Higgs portal, to have a vector portal interaction, 
we assume a U{iy gauge theory in the dark sector as the simplest model including 
a gauge boson. We also assume that only the SM quarks (but not the leptons) are 
charged under the U{1)'. In other words, we are dealing with a leptophobic vector 
portal interaction. 

The total Lagrangian consists of the standard model part, the dark sector and 
the interactions between these two sectors: 


C. = Tsm + ^DM + -^int ■ (1) 

The SM covariant derivative acting on the quarks must now be modified as 

( 2 ) 




where z is the dark charge of the quark field that the covariant derivative acts on. 

The dark matter Lagrangian consists of a fermionic dark matter and a complex 
scalar field both charged under C/(l)h 


/:dm = + X -m^)x+ , (3) 
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where is the 17(1)' held strength, x is a Dirac fermion as the dark matter 

candidate and 0 is a complex scalar held. The dark sector covariant derivative is 
given by 

D'^ = d,-ig'^Z'^. (4) 

We assume that the interaction of the standard model particles with the U{iy gauge 
boson, Z' is leptophobic. In other words, non of the leptons in the SM are charged 
under C/(l)h The Tint consisting of the scalar-Higgs and Z'-quark interactions reads: 

Tint = -A(#*) + s'^Z'^Ql^Ql + g'^Z'^UR^UR + g'^Z'^dR^dR , (5) 


where H is the SM Higgs doublet, Ql, ur and dR are respectively left-handed quark 
doublet, right-handed up-quark singlet and right-handed down-quark singlet. The 
couplings of the light quarks u,d,... with Z' are considered to be negligible. Therefore, 
by zq^ , Zuji and we mean the dark charge of only the third quark family i.e. the 
t and b quarks. The upshot is that the dark matter scattering off nuclei lacks the 
tree level contribution (see Fig. H and remains suppressed as expected from direct 
detection experiments. 

Having introduced dark gauge boson, Z', interacting with SM and DM fermionic 
currents, one should note that new anomalies from triangle Feynman diagrams may 
arise. However, it can be shown that assigning appropriate dark charges for quarks 
can lead to a anomaly-free theory 0. 

The anomaly-free conditions put constraints on the top and bottom quark 17(1)' 
charges: zq^ = —2, = +2 and = +2. Substituting these charges in Eq. ([5|) 

the Tint becomes: 


Tint — Ai 


(#*) [HH^)+g'Z'^irjH + g'Z'^brj^b. 


( 6 ) 


The Tint consists of the Higgs portal where the scalar field interacts with the 
SM Higgs quadratically, and the vector portal where the dark gauge boson interacts 
axially with the third family quarks. 

The scalar field (p does not interact directly with the DM particle y but interacts 
with that only through another mediator of the model i.e. the dark gauge boson 
Z'. On the other hand, the gauge boson mediator interacts directly with the dark 
matter particle which can be seen by following the red line in Fig. [TJ The scalar held 
(p and the gauge boson Z' are connected to the SM respectively via the Higgs portal 
and via an interaction with the third family of quarks (vector portal). The novelty 
of the current minimal model is that there are two distinct mediators at the same 
time which makes a bridge between the DM sector and the SM sector. Schematically 
these interactions are shown in Fig. [TJ 

^For some details on anomaly-free conditions in the extended SM including Uil)' interactions with no 
additional fermions see [46) . In our model we have an additional fermion that is the dark matter Dirac 
field. Taking equal dark charges for left- and right-handed components of the Dirac fermion we 

lead to the same anomaly-free conditions mentioned in [46) . 
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Figure 1: Two-portal SM-DM interactions: the dark matter candidate indirectly can 
interact with standard model through a Higgs and a vector portal. In this model the dark 
matter can interact directly with only one of the mediators i.e. the dark vector boson. 


The SM Higgs potential is given by 

Vh = -pn - Xh , (7) 

where the Higgs doublet takes on a non-zero vacuum expectation value (vev), 



0 

V + h 


We assume that the scalar mediator also takes a non-zero vev, 


{ 4 >) 


= v' + -^h'. 

V2 


( 8 ) 

(9) 


h and h' are respectively the SM Higgs and the singlet scalar field fluctuations around 
their vacuum expectation values. It is worth mentioning that once the complex scalar 
develops a non-zero vacuum expectation value, v\ this breaks the U{iy symmetry 
spontaneously and the Z' boson will acquire mass. On top of that, the strength of 
these vertices, Z'Z'h and Z'Z'h', is proportional to v'. Therefore, according to the 
Feynman diagram in Fig. [2l the choice = 0 will give rise to zero DM-quark elastic 
scattering cross section which is a trivial scenario. 

After substituting Eq. Q in Eq. ([3|) and expanding the Lagrangian, the mass 
of the dark gauge boson turns out to be g'v' 1^/2. As may be followed in [16] the 
masses of the SM Higgs particle h and the scalar mediator h' can be obtained by 
diagonalizing the mass matrix. 


M = 


h 

h’ 


h 


f 2 \hv^ 
V y/2\ivv' 


h' 

\/2Xivv' \ 
— ^Xiv^ ) ’ 


( 10 ) 
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where we have used the following relations coming from minimizing the total poten¬ 
tial, 

ml =- Xiv^ , ( 11 ) 

- Xiv'^. (12) 

We can redefine the scalars h and h' by introducing a mass mixing angle in order to 
get a diagonalized mass matrix, 

h = sin {6) h + cos [9) h! , (13) 


h' = cos (0) h — sin (0) h ', 


(14) 


with the mixing angle 0 being 


tan (0) 


1 

1 ’ 


V = 


mi — mi, 

h h' 


where is the off-diagonal entry of the mass matrix in Eq. m- 
The masses of the redefined scalar fields read, 


(15) 


2 2 


m? -b m? 
h h' 


mi — m 
h 


h' 




(16) 


where the upper sign (lower sign) corresponds to ruh {my). The standard model 
Higgs is denoted here by h with mass mh = 125 GeV and h' is the singlet scalar. 
Exploiting Eq. (fT^ and Eq. (fT6l) we can obtain the quartic couplings as a function 
of SM Higgs mass, singlet scalar mass, the mixing angle and the vacuum expectation 
values V and v', 


X 


Xh = 


m\i sin^ 0 -b cos^ i 


m\i cos^ 9 X-ml sin^ 0 


n2/2 
Ai = 




mf — mil 
^ ^ sin 20 , 


2'/2vv' 


(17) 


The vacuum stability of the total potential is equivalent to having positive eigenvalues 
for the scalar boson mass-squared matrix. At tree-level this brings in the following 
constraints on the couplings (see also [E]); Xh > 0, and v'‘^{XhX—2Xi) > 

v^XiXh- In our numerical analysis we will choose m^,my, 0, v' and g' as free 
parameters. 


6 











3 Invisible Higgs Decay 

In the present model, there are two new decay channels for the SM Higgs which 
can modify the total decay width of the Higgs boson within the SM. The current 
measurement of total decay width for the 125 GeV Higgs reads, ~ 4 MeV m- 

In case the dark gauge boson is light enough such that mz' < m,h/2 the Higgs boson 
is kinematically allowed to undergo the following invisible decay 

rinv(h Z'Z') = ^ ^ (1 - 

lOTTnih 

In addition when we consider light scalar boson with rrih/ < ruhj^, another decay 
channel is plausible for the SM Higgs with 

Tin^ih h'h') = (19) 

where 


c = 3V2Xv' cos^ 9 sin 6 + 12Xhv cos 9 sin^ 9 — 6Aiu cos 9 sin^ 9 (20) 

+2Aiu cos 9 + gV^Xiv' sin^ 9 — Ay/^Xiv' sin 9. 

We thus expect the total Higgs decay width to modify as 

rffiggs = cos^ 9 r|^gg + e{mh - 2mz')T{h ^ Z'Z') + Q{mh - 2mh')T{h ^ h'h'),{21) 

where 0 is the step function. It is worth mentioning that since here Higgs has no 
tree level interaction with the DM, the invisible Higgs decay width is independent of 
DM mass at leading order. There exist an experimental upper limit for the invisible 
branching ratio of the 125 GeV Higgs decay investigated at the LHC, Brjnv ^ 0.35 
|48j . In our numerical analysis when applicable, we restrict ourself into the parameter 
space which satisfies the condition Linv/LH^gg < 0.35. 


4 Direct Detection 


The tree level DM-quark elastic scattering is suppressed because dark vector boson 
interaction with light quarks are assumed to be negligible. As depicted in Fig. [2] the 
hrst leading contribution to the elastic scattering amplitude is obtained through a 
one-loop interaction coupled to the SM Higgs or the Higgs-like scalar where the DM 
particle and dark gauge boson run in the loop. The DM-quark scattering amplitude 
is obtained as 


A4 = —4ig"^v'-^ 


mn 


sm 


cos 


qq X 


V \-{pi - P2)‘^ - ml {pi - P2)‘^ - ml, 

_ x{P2)lii{^ + m^)7^x(4>i) _ 

(27r)2 [{P 2 - g)2 - mA][(pi - g)2 _ m|,][g2 - 


( 22 ) 
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Figure 2: The Feynman diagram for the DM elastic scattering with quarks. The wavy 
lines stand for the propagation of the Z' boson. 


where respectively pi and p 2 are the four-momenta of the incoming and outgoing DM, 
and the DM four-momentum in the loop is denoted by q. Since {pi —p 2 )^ ^ 
we can then perform the loop integral at t = (pi — ^ 2 )^ ~ 0 to get the effective 
scattering amplitud^ 


MeS = 


47r2mx V 


cos^ 6 sin^ 0 1 


m 


h' 


ml 


{qq){xx) = otq {qq)ixx) 


(23) 


where 


S{P) 


-2 + /31og/3 


/ 3 ^- 2 / 3 - 2 ^_ + 

x//32 - 4/3 


(24) 


and /3 = In order to find the DM-nucleon elastic scattering cross section 

one needs to evaluate the nucleonic matrix element. However, at the vanishing mo¬ 
mentum transfer we can make use of the conventional assumption that the nucleonic 
matrix element with quark current is proportional to the nucleonic matrix element 
with nucleon current |49H51j 


in which 


Y,aq{Nf\qqm = aN{Nf\NN\N,) , 
q 


OTV — 


E fr. 

q=u,d,s 




nin 


+ 


27 


^Tg Y1 




q=c,b,t 


rUn 


(25) 


(26) 


The scalar couplings and are responsible for the low energy strong interaction 
and nucleon mass is denoted by mw- In the numerical computation in Sec. [B] we shall 
use the following values for the scalar couplings [52] 

/P = 0.0153, /f = 0.0191, /P = 0.0447. (27) 


Spin-independent (SI) total cross section of DM-nucleon elastic scattering is finally 
achieved as 


^si —-1- 

TT 


(28) 


where is the reduced mass of the DM-nucleon system. 


^Currently, the typical recoil energy under examination at direct detection experiments is En ^ 10 
KeV. On the other hand, the momentum transferred to a nucleus of mass Mjv is given by t = 2MpfEp>. 
For a xenon nucleus for example, we obtain t ^ 2 x 10“^ GeV^. Thus we expect t <C 
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Figure 3: The Feynman diagrams for the DM annihilation processes. 

5 Relic Abundance 

The fermionic dark matter candidate in the model laid out earlier is of WIMP type 
DM whose present day density, the so called relic density, is a remnant from freeze-out 
epoch in the early Universe. The freeze-out mechanism is based on the assumption 
that dark particles had been in thermal equilibrium in the early time at temperatures 
T > rriDM- In an expanding Universe the annihilation rate of dark particles into SM 
particles slows down and there is an epoch with T <C ut-dm after which this rate 
descends below the Hubble expansion rate. On the other hand, from this time on 
dark particles are not kinematically allowed to get reproduced. Thus, in effect, the 
number density of dark particles, n^, remains asymptotically constant within the 
comoving volume. 

The leading DM annihilation reactions which are necessary to determine the 
relic density are shown in Fig. [3l In this figure the first and third annihilation 
processes occur via a Z' boson exchange in s-channel: XX bb,tt, Z'h, Z'h', while 
the second diagram shows annihilation with an intermediate DM via t- and u-channel: 
XX Z'Z'. The Boltzmann equation provides us with the evolution of DM number 
density in terms of thermal averaged annihilation cross sections (o'ann^’rei) as 

+ 3Hn^ + {(Ji,nnVre\)[n\ - = 0, (29) 

where, is the DM number density at equilibrium condition and H is the Hub¬ 
ble parameter. In order to determine the present value of the number density and 
therefore the relic density one should solve numerically the Boltzmann equation at 
freeze-out condition which is when the dark particles are away from equilibrium. 

We first implement our model into the program LanHEP [53] to give us all the 
basic vertices and Feynman rules of our model. Later on to analyze the DM relic den¬ 
sity we employ the package MicrOMEGAs |52j which requires our output files from 
the LanHEP program. To check the validity of our model implementation into Lan¬ 
HEP, we utilize the program CalcHEP |54j using our LanHEP outputs to calculate 
the annihilation cross sections. From this we found agreements with our analytical 
calculations given in appendix for the relevant annihilation cross sections. 


6 Numerical Analysis 

In this section we will find the viable region in the parameter space which respect 
observed relic density, invisible Higgs decay width measurement and constraints from 
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Figure 4: Results for the spin-independent elastic scattering cross section of DM with 
proton as a function of two variables, scalar vacuum expectation value v' and the mixing 
angle 9. The coupling g' is generated in the range 0.001 <g' < 1. For both values of the 
DM mass, the elastic scattering cross section is below the upper limit given by LUX and 
XENONIOO. 


direct detection experiments. We consider in our parameter space as independent 
free variables: m^, nih', g', v' and 9. Throughout our study we keep fixed the 
SM Higgs mass as ruh = 125 GeV and the SM Higgs vacuum expectation value as 
V = 246 GeV. As a first numerical look, we would like to find the viable region in the 
parameter space for a given set of variables by scanning over the scalar 

vacuum expectation value v' and the mixing angle 6. To do so, we pick out random 
values for the coupling g' in the reasonable range of 0.001 < < 7 ^ < 1. Our results 
presented in Fig. [4] for two sample values of the DM mass, indicate that it is possible 
to find viable regions in the parameter space for the wide range of 100 < v' < 1000 
and the mixing angle 0 < 0 < 7 r/ 2 . 

To move on, we fix two variables out of five independent free parameters when we 
scan over the parameter space. For the vacuum expectation value of the singlet scalar 
we choose v' = 800 GeV and perform our calculations for two different mixing angle 
with sin0 = 0.01 and sin0 = 0.1. It is then ensured that with these choices and the 
range of the masses we will pick out for the singlet scalar, the quartic couplings will 
respect bounds from perturbativity and vacuum stability conditions when relations 
in Eq. (HZI) are applied. 

We begin our scan over the parameter space by generating random values of 
order ~ 10® for three free parameters in the ranges: 1 GeV < < 1 TeV, 20 GeV 

< m/i/ < 150 GeV and 0.01 < < 7 ' < 1. Given the mass relation mz' = g'v'jy/2, 
the Z' boson mass will then lie in the range 5.6 GeV < Z' < 565 GeV. We then 
use the combined results from Planck [T] and WMAP [2], 0.1172 < Q.h? < 0.1226, 
to exclude large regions in the parameter space which are inconsistent with these 
observations. At the same time, when mz',mhi < mhl2 we check further to make 
sure each generated point in the parameter space can fulfill the upper limit constraint 
on the invisible Higgs decay width. 
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Figure 5: Spin-independent elastic scattering cross section of DM with proton is shown 
as a function of DM mass. All the points displayed in the plots respect constraints from 
observed relic density and invisible Higgs decay width measurement. The anticipated 
upper limit bounds on the elastic scattering cross section imposed by LUX and XENONIOO 
experiments are placed to make comparison. The mixing angle is chosen sin0 = 0.01 in 
the left panel and sin0 = 0.1 in the right panel. 

Using the formula provided by Eq. (p5|l we compute the DM-proton elastic scat¬ 
tering cross section in terms of DM mass within the parameter space restricted by the 
observed relic density and invisible Higgs decay width measurement. For a wide range 
of the DM mass our results for the elastic scattering cross section are summarized in 
Fig. [5] for sin0 = 0.01 and sin0 = 0.1. 

It is evident from these figures that in our model there exist viable regions in 
the parameter space with DM elastic scattering cross section well below LUX and 
XENONIOO bounds when the ratio mz'Imy < 5. With sin0 = 0.01, all the points 
with correct relic abundance in the left panel of Fig. [5] respect the upper bound from 
the invisible Higgs decay. However, for larger mixing angle sin0 = 0.1, the invisible 
Higgs decay constraint excludes some portion of the region with correct relic density 
as can be seen by the right panel in Fig. [3 

Therefore we emphasize here on an interesting feature of the two-portal model 
that the dark particle can evade direct detection in the range of DM mass from a few 
GeV up to 1 TeV. 

7 Gamma-Ray Emission Prom DM Annihila¬ 
tion 

The gamma-ray excess observed in the GC from the analyses of the Fermi Large Area 
Telescope (Fermi-LAT) data is one of the places to look for the trace of the dark 
matter signals. Among other disfavored scenarios such as millisecond pulsars and 
cosmic-rays sources, the annihilation of the dark matter (which is more accumulated 
in the center of the Galaxy) into SM particles explains well the observed gamma-ray 


11 















Figure 6 : The spectrum of the gamma-ray produced by DM annihilation for masses tudm = 
20, 34 and 44 GeV with slop parameter 7 = 1.15,1, 3 and 1.3 respectively for each DM mass. 
The DM mass tudm = 34 GeV is more compatible with Fermi-LAT data analysis [TO] . 

excess. 

After it was worked out in [5] where the excess reported for the first time, more 
accurate analyses were implemented by different groups confirming the original re¬ 
sults [Hi nni iHH] . 

In this section we examine the two-portal model discussed in the last sections 
for the gamma-ray excess. The region of interest (ROI) we use in our computation 
is that of considered in m, i-e. at Galactic latitudes 2° < |6|< 20° and Galactic 
longitudes |f| < 20° known as Inner Galaxy. 

Let us briefly review the material we use to obtain the gamma-ray spectrum from 
dark matter annihilation. 

The flux of the gamma-ray produced by annihilation of dark matter into SM 
particles is given by 




{av) dNy 

Svrmpjyf dE^ 


I 

J 1.0.S 


p^{r)dl, 


(30) 


where {av) is the velocity averaged total annihilation cross section, tu-dm denotes the 
mass of the dark matter and dN^ / dEr^ is the gamma energy spectrum produced per 
annihilation. The integral of the density squared is performed over the line-of-sight 
(l.o.s). The dark matter density as a function of r, the distance from the center of the 
Galaxy, is given by p{r). This density function is assumed to be spherically symmetric 
and is given by the generalized Navarro-Frenk-While (NFW) halo profile 


p{r) = po 


{r/vs 


1-7 


(1 -t- r/vsY ’ 


(31) 
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with the local dark matter density = 0.4 GeV/cm^ and the radius scale = 20 
kpc. Due to high uncertainty in the dark matter density near the center of our Milky 
Way galaxy, the inner slop parameter takes values in the range 7 = 1 — 1.3. 

We use micrOMEGAs to compute the gamma-ray spectrum for dark masses 
moM ~ 20,34,44 GeV that are picked out from the viable parameter space obtained 
in the previous section (see Fig. ED. The gamma slop parameter is chosen 7 = 1.15 for 
moM = 20 GeV and 7 = 1.3 for moM = 34,44 GeV. The annihilation cross section 
that we obtain for different masses are {av) = 2.14 x 10“^® cm^ s“^ for tudm = 20 
GeV, {av) = 2.41 x 10“^® cm^ s“^ for moM = 34 GeV and {av) = 2.34 x 10“^® 
cm^ s“^ for rriDM = 44 GeV. In Fig. El we have plotted the energy spectrum of the 
gamma production for the masses mentioned above. As seen in this figure the DM 
mass moM = 34 GeV has a better agreement with the Fermi-LAT data analysis 
performed in m- 


8 Conclusions 

In this paper we have proposed a minimal fermionic dark matter model with two 
Higgs and vector portals in tandem, both charged under a U{1)' symmetry. The Z' 
dark boson beside coupling to the fermionic DM has leptophobic interaction with the 
SM particles. The Z' coupling to the light quarks has also been considered negligible, 
resulting in a suppressed DM-nucleon interaction. The leading contribution to the 
DM-nucleon elastic scattering comes from a one-loop Feynman diagram as shown in 
Fig. El 

An interesting result is that we find a wide range of DM mass from a few GeV 
up to 1 TeV in the viable parameter space respecting all constraints from observed 
relic abundance, direct detection bounds and invisible Higgs decay width. 

The dark matter annihilation into SM particles results in the gamma-ray energy 
spectrum that fits best with the Fermi-LAT data for the DM mass ttidm = 34 GeV. 

The new aspects of the two-portal DM in comparison with the earlier works which 
involve either the Higgs portal or the Z' portal are discussed below. 

1) In the fermionic DM model introduced in [58], DM interacts with the SM par¬ 
ticles with exchanging the SM Higgs {h) or a new scalar particle {h'). They found 
out that almost the entire parameter space is excluded by the XENON and GDMS 
except a very small resonant regions corresponding to very restricted conditions for 
the DM mass: rriy, ~ mh/2 or rriy' ~ rrih 'It is likely that with future experiments 
like XENONIT, these small viable regions get even smaller. Moreover, the model 
discussed in [58] predicts a large invisible Higgs decay width, such that with the cur¬ 
rent bound from the LHG, the DM masses less than about 50 GeV are excluded. In 
contrast, the presented model in this paper with simultaneous two portals provides us 
with a rather wide viable parameter space well below XENONIOO and LUX bounds. 
More precisely, we see that the viable range for the DM mass is from a few GeV up 
to 1 TeV. It is quite unlikely that the future proposed or planned direct detection 
experiments can rule out the entire viable parameter space for the aforementioned 
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range of the DM mass. One more interesting feature of the proposed model is that 
the total Higgs decay width including the effects of the new particles does not de¬ 
pend on the DM mass. Therefore, even for small DM mass it becomes possible to 
hnd viable regions in the parameter space which respect the invisible Higgs decay 
width bound. 

2) In the two-portal model, the coupling between the fermionic DM and the 
Z' boson is of vector type while the anomaly free condition imposes a pure axial 
coupling between the Z' boson with the third family of quarks. One crucial feature 
of the present model is the fact that with the help of scalar mediator (Higgs portal) 
and the non-universality of the Z' boson, spin-independent elastic scattering of the 
DM with nuclei is possible. In the earlier works (see e.g. |26p27|) with the same axial 
coupling between Z' boson and the third family of quarks there is no DM-nucleon 
elastic scattering at all. Therefore we have shown in the two-portal model that even 
without having any Z' interaction with u and d quarks, there is a prospect for a 
signal at the direct detection experiments. 


A Dark Matter Annihilation Cross Sections 


We give the DM annihilation cross section formulas in this section. The annihilation 
cross section for the process XX ff with f = b,t is obtained as 


0-annt’rel(XX ^ //) = 


g' — 4mj /s _ 8777 , 2777,2 _|_ 2sm^ — ^sm'j- 


2Trs 


(s - m |,)2 m|,r|, 


(32) 


The next annihilation process XX hZ', is again mediated by a dark gauge boson 
via s-channel. We get the following result for the cross section as, 


0-annl^rel(XX ^ hZ') = 


g'%'^ sin2 0 (s + 2m2)^[l - (m^ + m|,)/s]2 - 4mlm%,/s^ 


IOtts (s — m\iY + 

We then obtain the DM annihilation cross section for the process xX h!Z', 


.(33) 


. 5'%'^ cos2 0 (s + - (K' + - 4m2,m|,/s2 

^ ^ - (,.-, 4 ,)^ + m|,r|, -■ 

Finally, we get the annihilation cross section for the process xX which takes 

place by mediating a DM via t- and u-channel. 


i^^rei(xx z Z) = 


87r2i 


dO. 


sm?^, — + \srri^ — 2m^ 

{t — — 777,2) 


(ttt,^ -|- ? 77 |, — t)^ + ts — sm^ + 2tm'^ -|- Am^rn?^:' + 2 ™-^ 

2{t — 7772 )2 


/OOnO O O 00 A 

ym't + m^i — u) -|- tts — sm't -t- 2um^ + + 2777^ 


2(77 — 7772)2 

where, s, t and u are the relevant mandelstam variables. 


( 35 ) 
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